The Ridge from the BFKL evolution and beyond 



Eugene Levin^'^''^ and Amir H. Rezaeian^ 
^ Departamento de Ftsica, Universidad Tecnica Federico Santa Maria, 
Avda. Espana 1680, Casilla 110-V, Valparaiso, Chile 
^Centra Cientifico-Tecnologico de Valparaiso, Casilla 110-V, Valparaiso, Chile 
Department of Particle Physics, Tel Aviv University, Tel Aviv 69978, Israel 
(Dated: January 21, 2013) 

We show that the long-range rapidity correlations between the produced charged-hadron pairs 
from two BFKL parton showers generate considerable azimuthal angle correlations. These correla- 
tions have no 1/Nc suppression. The effect of gluon saturation on these correlations are discussed 
' and we show that it is important. We show that a pronounced ridge- like structure emerges by going 

, from the BFKL to the saturation region. We show that the ridge structure at high-energy proton- 

proton and nucleus-nucleus collisions has the same origin and its main feature can be understood 
due to initial-state effects. Although the effects of final-state interactions in the latter case can be 
non-negligible. 
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I. INTRODUCTION 



O |. The main objective of this paper is to understand the long range rapidity correlations of charged-particle pairs in 
JL ' the azimuthal angle separation between the two particles around the near side Aip sa 0, the so-called ridge which has 
Qj', been recently observed at the LHC in -Js = 2.76 TeV Pb-|-Pb collisions [l| and also in ^/s ~ 7 TeV proton-proton (pp) 
(~| ■ collisions [3] . The CMS collaboration Q recently reported that the ridge type structure exists in pp collisions at ^/s = 7 
TeV for high multiplicity iV > 90 event selections. The origin of the ridge in pp collisions at the LHC is not still well 
understood and it has been a subject of growing interest, see for example Refs. [3|-l5|- The ridge was previously seen at 
_ RHIC in central Cu-|-Cu coUisions at y/s = 62 A GeV and in Au-|-Au collisions at y/s = 200 GeV 0]. The description 
ly-j , of nucleus-nucleus (AA) collisions is generally more complicated compared to the case of pp collisions. However, given 
[ — ' the relative similarity of the observed ridge structure in both pp and AA collisions in terms of multiplicity, transverse 
, momenta and rapidity separations of pairs, it is natural to ask whether the ridge phenomenon has a unique origin 
■ and can be understood only by initial-state effects. We recall that the highest multiplicity events per unit rapidity in 
' pp collisions at y/s = 7 TeV is compatible to the one in central Cu-|-Cu collisions at RHIC. 

In high density QCD, we expect large rapidity correlations for produced hadron pairs with the value of their 
transverse momenta about the gluon saturation scale Qs [1, HHi see also Ref. [9|. At first sight, these correlations 
should be small at fixed impact parameters. It has been argued [sj that in the color-glass-condensate (CGC) approach 
(lo| there is a source of the long range rapidity correlations which transforms into the azimuthal angle correlations 
due to the collective flow in the final state. In Ref. Q it is argued that such mechanism can qualitatively explain the 
azimuthal angle correlations in proton-proton collisions without a significant flow effect. The issue of the importance 
of final-state and collective flow effects in the observed ridge structure in pp collisions is still debatable [3] , see also 

Refs. [iil,li3- 

In this paper we will introduce a new source of long-range azimuthal correlations for the produced charged hadron 
pairs. We show that the intrinsic long-range rapidity correlations between the produced hadron pairs from two parton 
showers generate considerable azimuthal angle correlations which do not depend on the interaction in the final state, 
and because of this, these correlations have the same origin both in pp and AA interactions at high energy. These 
correlations have no l/Nc suppression as one considered in Refs. [1, Q. Recently Kovner and Lublinsky in a very nice 
paper fisj put forward a general discussion toward understanding the ridge. We have an additional goal here, like 
Ref. 13] we shall try in this paper to understand the general feature of the ridge based on very general grounds and 
will show that the main features of these correlations both in rapidity and emission angle can be simply understood 
within the BFKL Pomeron calculus [l3 - [l9| . The extension beyond this framework inside the saturation regime will 
be also discussed. 

The paper is organized as follows: In sec II, we introduce our mechanism for the azimuthal correlations and illustrate 
the main idea within the perturbative framework. In Sec. HI, we consider double inclusive gluon production and 
its correlations within the BFKL Pomeron approach. We show that the azimuthal correlations between produced 
hadron pairs from two BFKL parton showers have long-range nature and will survive the BFKL leading log-s re- 
summation. In Sec. IV, we provide estimates of azimuthal correlations in both pp and AA collisions in the BFKL 
and the saturation regions. As a conclusion, in Sec. V we highlight our main results. 



FIG. 1: Mueller diagrams for two parton showers production. The wave lines denote the BFKL Pomerons. This is the typical 
diagram which gives an angular coUimation about /S.tp ~ 0. 



II. THE AZIMUTHAL CORRELATIONS: THE ORIGIN 



In this section, we show that the long range rapidity correlations in azimuthal angle separation between the hadron 
pairs can be simply understood in the perturbative QCD approach. In the parton-like language, the Mueller diagram 
[2Qi l shown in Fig. [1] (right panel) describes the emission of two particles (partons) from two parton showers. One 
can write the contribution of this diagram to the cross-section of double inclusive gluon production in the following 
generic form, 

d<J 1 /" ,9 ,,9 x da / X \ da 



dyid^pi dy2d'^p2 



where Njph is the scattering amplitudes for Pomeron (ladder)-hadron productions along which transverse momentum 
Qt is transferred and da / dyid^PiT denotes the corresponding cross-section of the gluon production with rapidity yi 
and piT in each of the BFKL Pomeron ladders. This factorization is based on the leading Log-s approximation i gnoring 
enhanced Pomeron diagrams. Eq. ([!]) can be motivated [2l| - |25j using three main ingredients: Gribov Reggeon 
and Pomeron [14 18] calculus, AGK cutting rules 26] and Mueller generalized optical theorem '20']. In the Pomeron 
calculus the amplitude Nph is a new ingredient which can be written in the following form, 

Nph {Qt) = 9pn (Qt) + ^ 9ipp (Qt = 0) G^jp (Qt) {M^/s^) + (2) 

where n denotes the number of produced state with mass Af„ in the diffractive dissociation with M^ax as its maximum 
value (about 2 GeV) by which one can still express Np as a sum of resonances, gpn denotes the vertex of the Pomeron 
with this state {gpn — gpp for n = I) and Gsp denotes triple Pomeron vertex. The first term in Eq. ^ describe 
the contribution of the state with finite mass and this sum can be approximated by the sum of produced resonances. 
The second term is responsible for high mass contribution and can be described by the Pomeron contribution which 

leads to the factor (A'P/sq) ^ where AiP is the Pomeron intercept and so is the energy scale (so ~ IGeV) [2%, 
see Fig. [2] In the framework of the high energy Pomeron phenomenology it turns out that Qt dependence of the 
resonance contribution is much steeper than the one in the triple Pomeron term. In the BFKL Pomeron calculus 
this fact has a natural explanation: the resonance contributions are determined by the non-perturbative soft scale 
which is about 1 fm, while the triple BFKL Pomeron vertex has a natural scale of the order of the saturation scale 
which increases with energy. It should be stressed that Np^ has a very simple physical meaning, namely Np/^ is the 
probability to produce two parton-showers in hadron-hadron collisions. 

At first sight, one may expect that Fig. [T] describes two independent parton showers, and therefore there should not 
be any correlation between two produced gluons from these two parton showers. However, angular correlations stem 
from the Qt integration in Eq. ([T]). Due to this integration the contribution of diagram in Fig. [T]is not equal to the 
product of two single inclusive cross-sections leading to nonzero two particle correlation TZ ^ Q. In order to illustrate 
this simple fact, let us for the sake of argument assume that the gluon production cross-section in one parton shower 
is proportional to Qt • Pi.T, or in other words, 



(3) 




N(Qr) 





FIG. 2: Diagrams representing the Pomeron-hadron scattering amplitude Njph (Qt) as a sum of resonance contributions, triple- 
Pomeron diagram with the vertex denoted by Gsjp and etc, see the text for the details. The wave lines denote the Pomeron, 
while lines represent hadrons. 



In this case, Eq. ([T]) simply becomes 

dyid^Pi,T dy2d^P2,T J dyid^pi^x dy2d^p2,T ^ ' \ /V J 

= -Pi.T • P2.T {n/2) f dQ\ NU{Q\) {Q\) f , {Q\) . (4) 

J dyid^pi^T dy2d-^p2,T 

The above equation explicitly shows an angular correlation between two produced gluons in two parton showers. 
Having this equation in mind, in the next section we will explicitly show that the vertex emission of gluon from the 
BFKL Pomeron with Qt ^ (see Fig. ^ have a structure similar to Eq. (jH). 

For simplicity and clarity of the presentation, let us first work in the Born approximation, see Fig. [3]-a. In this 
approximation up to strong-coupling corrections, the inclusive singlet gluon production at very high energy, 
assuming that all components of the exchanged momentum are much smaller than the projectile and target momentum 
(for s >>| t I), is given by 



2a^CF f^,^ r,[qT,q'T)r''[-i<l-Q)T,-iq'-Q)T 



dyd'pT t:^ J'^'^^ ql{Q-q)'Tq'3{Q-^?T ' 

where Cp — {N^ — ^)/2Nc is the SU{Nc) Casimir operator in the fundamental representation with the number of 
color equals Nc- We used a notation pT — gr — q'-r- The effective vertex F^ and F^ for the emission of gluons (see 
Fig. |3]-a) are related to the Lipatov vertex Fj^^^ [iJ, in the following way, 

i^'' (qr^q'T) = ^Pi^Psi^rp^ (^(fT,g'T) > (6) 

where pi and p2 represent the momenta of the incoming projectile and target gluon, and the center of mass energy is 
s = 2pi.p2. The product of the two vertices appeared in Eq. ([5]) can be simplified to, 

K{QT;qT,q'T) = ^r^(?'T,gT) f^(-('f-0)T,-(9^-Q)T) , 

= 4(9t(Q-9)I + ql{Q-q'?T " pIQt)- (7) 
Pt ^ ' 

Substituting the above expression into the cross-section Eq. ([5]) , one immediately obtains 

oc I ^ + } ^ y\v\ ^ ^ \ 



dyd^PT J pI ^q'^{Q-q)l ql{Q-q% q^ q'3 {Q - q)l {Q - q') 



T 



> a 



2 4, ^ ^ 

Pt It It 



PT :>qT;QT <^qT^ ^3 f d (JT , ^Qt ' PT 



(pT ■ Qt^ 



{2 + 2^ + 4 V^:!^}. (10) 



2 4, ^ ^ 

qrp Prp I. Prp Pj, 



Notice that in the Born approximation we do not consider the kinematic region <Si Qt since we will show later 
that this region is not important for the azimuthal correlations from the BFKL Pomeron. Moreover, we should 
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Fig. Ob 



Fig. Ha 



FIG. 3: The ladder-type diagram that describes the production of gluon with transverse momentum piT in the Born approx- 
imation (Fig. [3]-a) and BFKL Pomeron (Fig. [3]-b) . The blobs represent Lipatov vertices and asterisks in left denote reggeized 
gluons. The produced gluon in the zth rung is shown within a box. 



stress that the expansion here are only for the purpose of illustration to trace back the origin of the azimuthal angle 
correlations in our approach while for the practical estimates, one has to perform the integrals without resorting to 
any approximation. 

First notice that Eq. ([5]) is symmetric^ under qt ^ <i't ^^'^ Pt ~^ ~Pt- In the expansion given in Eqs. (|9I10I) . 
we changed the variable to q'q^ = qr — Pt- Changing the variable in Eq. ^ to qr = q'r + Pt and then in the 
same fashion expanding we get the same expression as in the above equations but the second term in Eqs. ()9ll0p 
will be with the opposite sign. Actually, these two expansions correspond to different regions of integrand in Eq. ([S]) . 
Summing these two contributions^ we obtain the following form for the double inclusive cross-section from Eq. U) in 

the case of ^ qT,q'Ti Qt ^ q't^Qt, 

da /" .9 ^ . ,9 , „ . da , „ _ , da 



dyi dy2 d'^pi^T d'^P2,T 



d^QT Q\ Nl^ {Qt) , (Qt = 0) , (Qt = 0) , 

(11) 



where Ai^ denotes the angle between pi,t and p2.T and we defined 

da 



4—^— / d^qT(l){qT,-qT) <P{qT - PT.PT - qr) , (12) 
Cf Pt J 



and 



dyd^PT Cf Pt 



a'^^^ ^ f '^^A,(- - ^ - - - ^ / ^ da_ 

^— — 2 / -^'P{qT,-qT) (p{qT " PT,PT - qT) ^ (— ) ^..j2^ ■ (13) 



7 JO-* ^ ^ 2 1 At' \hi 1 H± J \H1 l'± 1 1' 1 ^J- J \ 4 / J J9 

dya-^PT Px J Qt Pt 

In the above, we used the fohowing notation, 

1 J -^(t)iqT,-qT) (kiOT -PT,PT - qr) 
/—\ = ^ 

^T j d?qT<t>{qT,-qT) <i){qT - pt.pt - qTY 

where denotes the unintegrated gluon density of the projectiles |27| for Qt ~ 0, 



(14) 



^ We thank our referee for drawing our attention to this point. 

^ These two expansions can be also envisaged as two different processes: in Eqs. Il9lf OH the transverse momentum of produced gluon is 
compensated by the gluon with the value of the rapidity smaller than the rapidity of the produced gluon with the transverse momentum 
Pt (gluon with rapidity in Fig. |3]-a), while expansion in q'j, we consider the process where pT is balanced by the gluon with the rapidity 
larger than the rapidity of the produced gluon with the transverse momentum pT (gluon with rapidity Y in Fig.[3}a). 
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(j){qT,-qT) = — 2- (15) 

Notice that at Q = the inclusive cross-section given in Eq. ^ is identical to Eq. (IT^ . 

This simple example indicates that we have a natural mechanism for the azimuthal correlations in the framework 
of perturbative QCD which does not depend on the final state interactions and leads to the correlations inside of 
initial wave-function of the incoming hadrons. In the next section, we will show that this azimuthal correlation has 
long-range nature and will survive the BFKL leading log-s resummation. 



III. LONG-RANGE AZIMUTHAL CORRELATIONS FOR TWO BFKL PARTON SHOWERS 

The generalization of the Born approximation to the case of gluon emissions from the BFKL Pomeron cannot be 
simply obtained via Eq. (|lip by replacing the unintegrated gluon density to the one obtained from the BFKL 
equation. Indeed, the unintegrated gluon density (j) depends also on Qt and we have to be very careful with putting 
Qt ~ 0. The inclusive gluon product can be generally written as 

rf^d^ " j lTK[QT;qT,q'T) ^,^^^_^^^ ^[Y-y,qT,QT-qT) ^ [y,qT - PT,Qt - qr + Pt) , 

(16) 

where K(^QT',qT,q'Tj the BFKL kernel given in Eq. ([7]) and we defined q'j^ — qr — Pt- In the above, the 

variable Y — ln(s/m^) denotes the total rapidity in the lab frame where m is the nucleon mass and y and px are 
the transverse momentum and rapidity of the produced gluon, respectively. Notice that at Qt = 0, the above 
expression has the same functional form as the fcy factorization (27j . The only dependence on pT comes from the 

term ^gV — Pt, Qt — qr + Pt^ for which we have the color-singlet BFKL equation (13. [28j: 

(t>{y,q'T,QT - qr) = — f f (fq" K(QT;q'T,Q"T] „2,^ — ^ <t) (vi+i, q"T, Qt - q"T 

^ J J Qrp {(J Q )t-- 



iff {Q - ?')!• 

where we defined as — aNc/TT. We first substitute (j) (y, q't^ Qt — 9t) given in Eq. pT)) into Eq. and expand the 
kernels of both equations up to the terms of the order of Q^. Then we again use Eq. ([T7| but at Qt = and collect 
all terms into 4> (^y, q'T, (Tr) ■ Therefore, we obtain the following equation. 



1^ - .If^/^..(0;f.,,v)^*(VW.-«H-V-,V){l.^..(^ 

+ • • • + terms of the order of Qt that do not lead to azimuthal angle correlationsj. (18) 

In order to understand better if the above approximation can be justified, let us examine the ladder summations 
which leads to the BFKL equation. At leading log-s approximation, the imaginary amplitude A of the quark-quark 
elastic scattering with exchange of a color-singlet gluon ladder whose vertical lines are reggeized gluons [13, [la 
can be written as 



^ ^ J qf+i.AQi+hT ~ QtY \p^+iJ 



glK[QT,qi,T,qi+i,T] f \ <^<3(9>+i,T)-(-eG(«,+i,T-QT) 

(19) 



where k{Qt-, q, <z'^ is again the BFKL kernel given in Eq. ([T7|) . The right-hand side of the above equation shows that 

the BFKL Pomeron can be written as a sum of production cross-sections as it follows from the optical theorem. The 
symbol denotes the integrations over n + 2-body phase space and the parameters f3i (with /3o = 1) is the standard 
Sudakov variables for the momentum of the i-channel gluons which obeys strong ordering of the longitudinal momenta 
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[ij, |15|, |28| . The expression in Eq. ([T^ takes into account the reggeization of gluons in t-channel that means that the 
spin of the gluon is not equal to 1 as in perturbative calculations but it is given by the reggeized gluon trajectory 

aa {q^.T) = 1 + eaiq^^T) = I + - I J'^'^'^h . (20) 

We recall that the produced gluon in the ith rung ladder is on-shell with pi,T = qi+i,T — qi,T- Then, in order to find 
Qt and Pi^T correlations one needs only to keep Qt ^ in the ith rung of the ladder (see Fig. [3]-b ) and to put 
Qt = in all other rungs. The contribution of this particular sell to the amplitude has the following structure 



K(^QT,qi,T,qi+i,T 



iqi+i,T - PiY{qi+i.T -Pi- QrYq^+i^Tiqi+i.T - QtY 

^. \ CG(<Ji+l,T)+eG(9i + l,T-QT) /p ^ \^ f-G{qi + l,T-Pi) + <iG{qi+l,T-Pi-QT) 



(21) 



Although the above equation includes the virtual radiative corrections, but has a very similar structure to the case 
of the Born approximation given in Eq. ([5]) and consequently in the same fashion discussed in the previous section, 
it also gives rise to the azimuthal correlations. Therefore, in order to extract the correlations between two produced 

gluons, it is sufficient to use Eq. ([T7|) in which we can put Qt = in (^q"T^ Qt — ^''t) ^^'^ 4' (qr^ Qt ~ q'rj ■ Using 
Eq. (jl8p and adding the contribution of the integration region in qx where {qr ~ Pt\ 3> \pt\, we obtain from Eq. ([T}, 

da f ,^9 ,^9> da ,^ „, da 



dyi dy2 d^pi^r <Pp2,T 



[ dQl NlMl) , (Qt = 0) , (Qt = 0) 

J dyi d^pi^T dy2 a^P2,T 

{l + IpItpItQt{^)' (2 + cos(2A^))}, (22) 
-^(l + IpItpIt ((Q^))(^)'(2 + cos(2A^))), (23) 



q 

where Aip is the angle between pi.T and P2.T and we defined the following notations, 



1 / </»(>"- y,qT,-qT) (l>{y,qT -PT,PT - qr) 

/ d^qripiY - y,qT,~qT) <P{y,qT - Pt.Pt - qr) ' 



(24) 



with n = 1,2. The normalization factor Af in Eq. ([23| is given by 

M ^ ^ fdQl Nl^{Ql) j^, [Qt = 0) , (Qt = 0) . (26) 

J dyi d^pi.T dy2 d^p2,T 

From the above, it is obvious that the production of two parton showers with a transverse momentum Qt along the 
Pomeron ladder, naturally leads to the long range rapidity correlation in azimuthal angle while the emissions from 
one parton shower given by the BFKL Pomeron contribution does not lead to such correlations, see also Ref. 0- 



IV. ESTIMATES OF AZIMUTHAL ANGLE CORRELATIONS IN pp AND AA COLLISIONS 

We recall that the long-range azimuthal angle correlations obtained by Eq. ([25)1 is valid in the leading log-s ap- 
proximation at high-energy. The azimuthal angle correlations in Eq. (|23p is uniquely determined by only knowing the 
average values (l/g|,") and ((Q|?*)). This equation was truncated at n = 2 assuming that the transverse momentum 
Qt in the Pomeron ladder is small. Let us explore the idea that Eq. is also valid in the saturation region (or at 
least on the boundary between the BFKL and the saturation regime) by choosing the corresponding average values 
{l/q^) and ((Q|")) in that region. 

In the kinematic regime of the BFKL (ignoring the saturation effect) from Eq. ([M]) we obtain « 
l/max{/i^", Q^"} where /i is the non-perturbative soft scale. At the LHC energies, the inclusive production stems 
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from the kinematic region in which saturation effects are important [2 91 - 13 1| . In this region, the interaction between 
Pomerons leads to more comphcated diagrams, the so-called enhanced diagrams shown in Fig. S) It has been shown 
(see Ref. [l^ and references therein) that the enhanced diagram leads to the value of the characteristic momentum 
of the order of Qs, namely we have Qt oc Qs ■ We do not need to follow complete calculations of this paper to 
understand why it happens so. Indeed, assuming that Qt ^ typical q « Qs we can replace the BFKL Pomerons in 
the loop by the Pomerons at Qt = 0. Therefore, in this case, we have (PQt = Qs- ^'^^ Qt ^ 9 the Pomeron 
exchange falls down with Qt making the integral being concentrated at Qt = Ht — Qs- In order words, if densities of 
partons in one parton shower is so large that we have already reached the saturation region of the gluon density, we 
can assume that the average (l/q^") » 1/Qs" where Qs is the saturation scale. This also follows from the high-density 
QCD within the CGC approach |1C)I | which describes the LHC data for the inclusive hadron production both in pp 
and AA collisions [2^, |30| , see also Ref. "sT] . We also assume that the density of partons in both parton showers is 
very large and consequently Pomeron enhanced diagrams are important, and therefore we can have ((Q^")) ~ Qs"- 
Therefore, we assume that in the saturation region we have only one relevant scale, the saturation scale, and the 
average transverse momenta are related to this scale. 

Notice that the maximum of the double inclusive production reaches at Pi,t ~ P2,t ~ Qs- Admittedly, we do 
not have a rigorous proof of this at our disposal without invoking any approximation, but this may be immediately 
understood within the CGC approach since Qs is the only dimensional parameter of the approach. This can be also 
seen in the simple case of the Born approximation by comparing Eqs. ((9|fT0|). Note that Eq. ([9]) gives the contribution 
at small values of pT and the correlations vanish at — >■ and increase with pT while Eq. ()10|) shows that the 
correlations falls down at large values of pt- Therefore, the correlation function has a maximum at pT ~ (qt), where 
{qt) is the typical transverse momentum of the system. The same argument is valid for the general case of the gluon 
pairs production from the BFKL Pomeron. This can be seen by comparing Eq. (j23p and its corresponding equation in 
the limit of pt ^ Qt ^ Qt- It should be stressed that the experimental data from the CMS collaboration indicates 
that the maximum of correlations occurs at the kinematic region that the saturation effects is important . 

The probability for the events with multiplicity equals N = 2{N) where {N) is the multiplicity in one parton 
shower, can be obtained by Eq. (1^51) and the corresponding cross-section of such events is (j{N — 2{N)) cx Af. Using 
Eq. (j23p . we obtain two-particle correlation function TZ for the event selections with multiplicity N as, 

7^(A^; 2/1,^2) = ''f^~''^'''%T 1 = + + cos(2A^))} - 1, (27) 

where the parameter n = E {N/{N)) is the relative average number of Pomeron parton showers in the event selections 

with multiplicity N, the average multiplicity (N) denotes the multiplicity in the mini-bias and function E gives the 

integer value of its argument. The pre-factor in Eq. (j27p comes from the counting the various possible ways to have 

two gluons production out of rl Pomeron parton showers. In other words, for simplicity we assumed that n showers are 

produced and two correlated gluons comes from only two different parton showers. The number of these pairs is equal 

-/- i\ 1 , d^N _ d^Af (one parton shower) , ^ i r • 

to n{n — 1) /2 and moreover we have ^y^i^. ^ = n dy d'^p- t ^' therefore the pre-tactor m Eq. (|27p can 

be readily obtained. 

Notice that the main background for the double inclusive gluons production is due to two jets production from 
one parton shower. However, this production is suppressed by making selection in the events. From AGK cutting 
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rules [23, it follows that the multiplicity in one parton shower is equal to the average multiplicity measured by the 
experiment in the mini-bias events. It should be stressed that the AGK cutting rules also work for two parton showers 
production in QCD . It is well-known that the gluon distribution in the BFKL Pomeron is close to the Poisson 
distribution, see Ref. [36[ and references therein. The production from two parton showers starts to be significant 
only for the events with multiplicity larger than 2(A^) where {N) is the mean multiplicity, see Fig. [T] On the other 
hand, the probability to have events with multiplicity 2(iV) in one parton shower is approximately suppressed as 
exp (-(2(iV) - {N)Y//2{N)) < 1 for the Poisson distribution. 

One can observe in Eq. ([27]) that except the over-all pre-factor, the coefficients does not depend on multiplicity 
and rapidity of pairs. Of course this feature may be altered due to possible contamination of two giuons production 
from one parton shower which may lead to short range rapidity correlations in the azimuthal angle A(^. However, 
in particular experimental set up with high multiplicity events where our underlying saturation assumption namely 
((*9t")) ~ ('Zt") ~ Qs" is at work, these correlations could be ignored and can only create a back ground that will 
fall off at large multiplicity events. 

In order to understand how much the azimuthal asymmetry depends on the value of ((Q^")), we next estimate 
((Qt")) in ths BFKL kinematic region ignoring the so-called enhanced diagrams (shown in Fig. 2]). In order to 
calculate ((Q^")) defined in Eq. (pS)) . we should know the non-perturbative amplitude Npp [Qt) defined in Eq. ©. 



For Npp (Qt), we use the quasi-eikonal approximation |24| . In this approximation we restrict ourselves to the first 



term in Eq. ^ and the contribution of the other terms is taken into account by introducing an extra factor Nq, 

NppIQt) = Nog^piQT), (28) 

where gpp is the vertex of Pomeron-proton interaction. This approximation has been widely employed in Pomeron 
phenomenology and works quite well in the description of the experimental data [s^ . The dependence of the BFKL 
Pomeron on the transverse momentum Qt is given by gpp (Qt) = l/(l-l-Q|i/m^)^ with the typical mass m determined 
from the experimental data. The dipole form of gpp is inspired by the Qt dependence of the electromagnetic form 
factor of the proton. Using this distribution we obtain (Qj) = m^/6 and (Qt) = The experimental data 

for diffractive production of the vector meson in the DIS indicates that = 0.8 GeV^ However, the CDF 



data on double jet production [3J| shows that the typical value of Qt could be larger leading to a bigger value 
for TO^ = 1.6 GeV^. A gain assuming that {qj^} = Q^" the corresponding two-particle correlation function TZ for 
Pit — P2T — Qs and n > 2 becomes, 

7e(A^;yi,y2) = ^^^^^^ { 1 + ^ (2 + cos (2A(p)) } - 1. (29) 

It is seen from above that the coefficients in TZ now depends on the rapidity via the saturation scale Qs in contrast to 
Eq. (P7)) . However, one should note that deep inside the saturation region the above equation is not reliable and one 
should then use Eq. (P7|) . It is instructive to notice that the two BFKL parton showers contribution lead to Eq. 
with the soft scale ^ instead of Qs- This scale is a new phenomenological parameter which does not depend on energy, 
and it is certainly fi < Qs- 

In Fig. [5] (right) we show the azimuthal correlation TZ obtained from Eq. ([29]) when {{Q'^)) was calculated within 
the BFKL region for two different masses m = 0.8 GeV^ and = 1.6 GeV^. In this plot, we take a fixed saturation 
scale Q^s — 0.6 GeV^. The chosen saturation scale is in accordance with the estimates of Ref. [2^ in pp collisions 
at the LHC. In Fig. [5] (right) we also show the azimuthal correlation TZ obtained from Eq. (|27p in the saturation 
region at different multiplicity N — n{N). It is observed that deep inside the saturation region we have the ridge-type 
structure, namely a second local maximum near Aip « independent of rapidity when pT is about the saturation 
scale. By comparing the results shown in Fig. [5] from Eqs. (|27I29|) . it is notably seen that by going from the BFKL to 
the saturation region, a pronounced ridge- type structure emerges. In Fig. [S] we show the experimental data from the 
CMS collaboration z| for projections of two dimensional correlation functions onto A(p (denoted in Fig. [S] left panel 
by A(j)) for the difference in pseudorapidity of pair 2 < A77 < 4.8 in different pt and multiplicity bins at 7 TeV pp 
collisions and reconstructed PYTHIAS simulations [13] ■ It is important to note that PYTHIA8 qualitatively fails to 
reproduce the local maximum in the near-side correlation in any of the pT or multiplicity bins [5], see Fig.E] Notice 
that our definition of two-particle correlation TZ defined in Eq. (|27p is different from the experimental definition R ^ 
shown in Fig. [5] with a over-all factor. Here given the simplicity of our approach we do not wish to compare directly 
our results with the experimental data. A meaningful comparison requires inclusion of correlations effect within one 
parton-shower, fragmentation and possible short-range correlation effects. Nevertheless, it is seen that the general 
feature of the near-side two-point correlations obtained by Eqs. I|27I29[) is compatible with the CMS experimental 
data m . One should note that for a denser system, the imposed condition of piT = P2T = Qs in Eq. (P7| shifts the 
relevant kinematic windows of the angular correlations to the higher pT since the saturation scale will be larger for a 
denser system. 
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A(p 

FIG. 5: Right: The correlation function TZ at diflerent multiphcity A'^ = n{N). The curves labeled by "Saturation" are the 
results from Eq. (|27|l when {{Qr^)) ~ Qs" with Qs being the saturation scale. The curves labeled by "Semi-saturation" 
are the results from Eq. (|29|l when ((Qt")) was calculated within BFKL region for two different masses m? — 0.8 GeV^ and 
m = 1.6 GeV^ In both cases we assumed (I/qt") ~ ^/Qs"- Left; Experimental data from the CMS collaboration for 
projections of 2-D correlation functions onto Atf) for 2 < Ar/ < 4.8 in different pr and multiplicity bins at 7 TeV pp collisions 
and reconstructed PYTHIA8 simulations [37| . Error bars are smaller than the symbols. The plot in the left panel is taken 
from Ref. 2]. 



Next, we consider the long-range correlations in nucleus-nucleus scatterings. It is straightforward to generalize 
Eq. p2p for the case of nucleus-nucleus collisions in the framework of the Glauber approach, namely assuming that 
multiple scatterings are only permitted on different nucleons while the nucleon-nucleon scattering stems from the 
BFKL Pomeron exchange. The double inclusive cross-section at a fixed impact-parameter between two center of 
nuclei b will then have the same form as Eq. except the extra dependence on the nuclear profile. Notice that the 
impact-parameter b is the conjugate variable to transverse momentum Qt- In the Glauber approximation for nuclei, 
the scattering amplitude of Pomeron-nucleus NpA in the region of small diffractive masses is defined as 

NpAiQr) = l^j cfbcfVe'"^^-^ SA{b-9)gp.p{9) 

- 9'pp{qT = 0) S\{Qt), (30) 

where gpp denotes the Pomeron-proton vertex and Sa (b) is the nuclear density profile defined by the Wood-Saxon 
parametrization. The second equation above is valid when the nuclear radius is larger compared to the proton size 
Ra >> Rp. Using Eq. ((30|) one can obtain the following expression for the double inclusive cross-section at fixed b 
in the framework of the Glauber approach in which the proton-proton scatterings are taken into account from the 
BFKL Pomeron: 

da A 1 daN duN 

dyi dy2 d^pi^T d'^P2,T d'^b 2 dyi d'^pi^T dy2 d^P2,T 

X {TlA{b)+ I pIt pIt {{^/qXrotonY (vlvlTlAib)) (2 + cos(2A^))}, (31) 
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where ^y^^-i^_^ ^ is the inclusive cross-section for proton-proton scatterings ^ and Taa is the nuclear overlap function 
for AA collisions. Using Eq. (|3ip one can calculate the correlation function TZ defined in Eq. ((27)) : 

7^(6,A^,2/l,y2) = ip?,TPlT((l/9')proton)'(v6V'7^lAW)(2 + cos(2A^)). (32) 

It is straightforward to show that in the Glauber approach, the inclusive production in AA collisions is proportional to 
the overlap function TAAib) while in the case of the double inclusive production instead it is proportional to T^j^{b). 
It is seen from Eq. ([5^ that independent productions are canceled in TZ at fixed b while in the integral over b the first 
term in Eq. pip gives the main contribution. It is worth mentioning that we do not need additional factor n as in 
Eq. (I29p since in the Glauber formulation for nucleus-nucleus scatterings, events with a fixed multiplicity correspond 
to a definite value of the impact-parameter. 

Deep inside the saturation region, the correlation function has the same form as for hadron-hadron collisions given 
by Eq. (P7)) with the saturation momentum replaced by that of the nucleus Ql{AA; x) ~ Taa {b) {pp; x) (3fl, |3l|. 
This is in agreement with the main idea of the CGC approach that difference between different reactions is only due 
to the different value of the saturation scale Qs- Assuming that we have (g^") proton = Qs" {pp',x) (n = 1,2) in the 
saturation region, it is seen from Eq. p2p that two-particle correlation TZ reduces by increasing the saturation scale of 
proton. However, the Glauber approximation is not reliable deep inside the saturation region, and one should instead 
use Eq. (P?]). consequently the slope of the reduction of the azimuthal correlations will be then different. One of the 
attractive feature of nucleus-nucleus collisions is that by using centrality cuts, one can study the underlying dynamics 
of two-particle correlations. 



V. CONCLUSIONS 



In this paper, we suggested a new mechanism for the long-range rapidity correlations in the azimuthal angle of 
produced hadron pairs, namely the long-range angle correlations of two parton showers component of the initial 
partonic (gluonic) wave-function. This mechanism can be conceived as a realization of the general ideas proposed in 
Refs. [1, [l3|. Our approach predicts large and of the same order long-range angular correlations both for hadron- 
hadron and nucleus-nucleus collisions inside the gluon saturation region. In our approach, the coUimation in Aip 
exists independently of the effects from flow in the later stages of the collisions. We showed that for extremely dense 
systems at the truncation level upto n = 2 for ((Qf-")) we have TZ — >■ 0.25 {TZ — —0.5 without correlation) at Aip w 0, 
and TZ still has a second local maximum near Aip « at ~ Qs- We showed that our mechanism qualitatively 
describes the main features of the observed ridge structure in proton-proton collisions at the LHC at y/s = 7 TeV. A 
detailed comparison with experimental data and numerical analysis is left for future. 

The main difference between our approach and the description in the framework of the CGC 043 is that in our 
approach the saturation region is explored from outside on the boundary with the BFKL region. We showed that 
a clear signal of the ridge-type structure emerges by going from the BFKL to the saturation regime. This is fully 
consistent with the fact that the saturation/CGC approach provides an adequate description of other 7 TeV data in pp 
collisions including the inclusive charged- hadron transverse- momentum and multiplicity distribution [29l Isoj . Finally 
notice that the correlations obtained in our approach is not suppressed with 1 /Nc in contrast to the prescription of 
Refs. d, 111 and survive in the leading order in 1/Nc expansion. 
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